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D , Abstract 



In Ref. (TJ Groote and Pivovarov have given notice of a possible fault in the use of sum 
rules involving two-point correlation functions to extract information on heavy quark pa- 
rameters, due to the presence of massless contributions that invalidate the construction of 
moments of the spectral densities. Here we show how to circumvent this problem through 
a new definition of the moments, providing an infrared safe and consistent procedure. 
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1 Introduction 

The vacuum polarization function suitable for extracting fundamental information of heavy 
quark-antiquark systems is built from the electromagnetic current j^(x) = cq Q(x) ^^ Q{x) of 
the heavy quark Q of mass M and eg electric charge : 

IV(g) = ijd 4 xe^ (0\Tj,(x)jt(0) |0> = (-g, u q 2 + q,q u ) I% 2 ) . (1) 

As it is well known two-point functions are analytic except for singularities at simple poles or 
branch cuts, the latter being originated by normal thresholds of production of internal on-shell 
states. Implicitly assuming that the absorptive part of II(g 2 ) starts at the massive two-particle 
threshold q 2 = AM 2 , vanishing below this point, the correlator satisfies the once-subtracted 
dispersion relation Q : 

n(g 2 ) = i% 2 ) - n(o) = q - — - } >- . (2) 
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This dispersion relation has been extensively used to determine heavy quark parameters within 
the method of sum rules because it allows to relate experimental input, on the right-hand side, 
with theoretical perturbative evaluations on the left-hand side 0. Indeed Imll(g 2 ) refers to 
the total cross section of heavy quark production a(e + e~ —>■ QQ). In practice, the spectral 
density Imll(s) is poorly known experimentally at very high energies and, in addition, we are 
interested in the very low energy domain because it is more sensitive to the heavy quark mass. 
Therefore one uses derivatives of the vacuum polarization at the origin, called moments, to be 
responsive to the threshold region : 




Mn = -A-rt) n(g% 2=0 . (3) 



Until present the evaluation of the perturbative two-point correlation function II(g 2 ) has 
only been carried out completely, with massive quarks, up to 0(a 2 ) [[§] and the procedure above 
has been termed consistent and effective in its task because the first branch point is set at the 
massive two-particle threshold. However in Ref. |l| Groote and Pivovarov have pointed out 
that at 0(a^) there is a contribution to the correlator which contains a three-gluon massless 
intermediate state (see Fig. [l|(a)). Its absorptive part starts at zero energy and, therefore, 
Eq. (^D is no longer correct. Moreover those authors have also warned about the fact that, at 
this perturbative order, the massless intermediate state invalidates the definition of the moments 
Ai n for n > 4 because they become singular. In Ref. an infrared safe redefinition of the 
moments, to cure the latter problem, has been provided; it consists in evaluating the moments 
at an Euclidean point q 2 = —se, se > 0, thus avoiding the singular behaviour. Nevertheless 
the fault in Eq. (^) due to the massless threshold still represents a problem because even if, as 
we will justify later on, we substitute this dispersion relation by 
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n(g2) = -_ _ <£t±± + 2_ _ 7^ , (4) 



1 Somctimes de Adler function defined as dH(q 2 )/dlnq 2 , to get rid of the subtraction constant, is used. 



a) 





b) 



Figure 1: (a) O(a^) diagram contributing to the vacuum polarization function of the heavy quark 
current (the vertical dashed line indicates the massless cut), (b) "Effective" diagram obtained by 
integrating out the fermion loops. It also has the topological structure of the "reduced" diagram 
that determines the massless cut singularity. 

(where the notation is self-explicative), the spectral function Imll3 ff (s) associated to the cut in 
Fig. 0(a) would hardly be implemented phenomenologically as gluons hadronize to both heavy 
and light quark pairs. Perturbatively the three-gluon cut would contribute to QQ production, 
i.e. to Imngg(g 2 ), but at higher order in a s . Therefore if we attach to an 0(af) sum rule 
analysis, that contribution should be extracted from the perturbative IT(g 2 ) evaluation. In this 
note we provide a bypass to recover the balance between the right-hand and left-hand parts of 
Eq. ©. 

2 Moments and the massless cut 



The perturbative contribution given by the diagram in Fig. |T](a) has been calculated at 
small q 2 (q 2 <C M 2 ) in Ref. 0. In this limit the quark triangle loop can be integrated out and 
it ends up in the diagram in Fig. [TJ(b) generated by an induced effective current describing the 
interaction of the vector current with three gluons 0, 
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dabc G*t G 
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G Pv 



(5) 



(6) 



where G% u is the gluon strength field tensor. The effective current in the QED case (G% v — > 
F^ u , a s — > a em , d a b c — > 1) can be easily identified from the lowest order Euler-Heisenberg 
Lagrangian (see Ref. 0). 

2 The permutations of the three gluons in Fig. fjl(a) are already included in the definition of the effective 
current. 



The correlator of the induced current (El) is then evaluated in the configuration space giving : 



(0\TJ,(x) Jt(0)|0) 
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(7) 



In momentum space we need to perform the Fourier transform of Eq. (|7|). Following the 
differential regularization procedure |J, which works directly in configuration space, the result 
for the vacuum polarization contribution of the diagram in Fig. |l](b) at small q 2 reads 
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with ft the renormalization point in this scheme, and d a b c d a b c = 40/3. 

As noticed by Groote and Pivovarov Jl|], moments associated to the diagram in Fig. |l|(b) 
are not defined if n > 4. Indeed differentiating Eq. (H) four times, at q 2 ~ 0, we get: 
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(9) 



whose real part clearly diverges if we set q 2 = 0. Larger n moments are also infrared divergent, 
and so the authors of Ref. (1| conclude that the standard sum rule analysis must limit the 
accuracy of theoretical calculations for the n > 4 moments to the 0(a 2 ) order of perturbation 
theory. 

One obvious way out of this infrared problem is to avoid the q 2 = point. As it has 
been discussed in Ref. |j], this solution is rather ill-conditioned from the phenomenological 
side though. Moreover we notice (as commented earlier) that it is not possible to implement, 
from the available experimental information, the second term in the right-hand side of Eq. (|J). 
However, if one does not insist in using full vacuum polarization for the sum rule analysis there 
is a way to overcome this infrared problem. 

3 Infrared safe definition of the moments 



The study of analytic properties of perturbation theory amplitudes shows that their singu- 
larities are isolated and, therefore, we can discuss each singularity of a perturbative amplitude 
by itself 0. As a consequence, II(g 2 ) in Eq. (0) satisfies a dispersion relation from Cauchy's 
theorem Q : 

„ 2-ni Js n s s 
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(10) 



Here [n(s)] n provides the discontinuity across a branch cut starting at the branch point s n . 

3 This expression also gives the residue Ri of a pole at s = Si if we interpret the discontinuity as ImIIi = 
irRi5(s — Si). However we do not consider the existence of QQ Coulomb bound states, as it is not relevant for 
our discussion. 



In the perturbative calculation, every discontinuity function [II (s)] can be associated to a 
"reduced" Feynman diagram obtained by contracting internal off-shell propagators to a point 
and leaving internal on-shell lines untouched. Its contribution is written down following the 
Cutkosky rules for the graph. The reduced diagram corresponding to the massless cut in 
Fig. 0(a) has the topological structure of the part (b) of that Figure. Let us emphasize though 
that our following discussion is not grounded on the q 2 <C M 2 regime where the fermion loops 
have been integrated out : the reduced diagram is just a symbol that specifies a singularity, 
and the black dots in Fig. 0(b) keep all the analytical structure of the fermion loops. 

In a general diagram the discontinuity across a specified cut needs not to be a pure real 
function in the physical region, only the sum of all cuts in a given channel gives the total 
imaginary part. Hence the separation between the imaginary parts coming from different final 
states, as performed in Eq. (f|) for the vacuum polarization, does not seem to come directly 
from the Cutkosky rules. Nevertheless in the heavy quark correlator the discontinuity across 
the three-gluon cut gives a contribution to the spectral function that is unequivocally real : 
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J dR 3g (0\f\3g)(3g\jl\0) , (11) 



from which the dispersive part can be evaluated independently of the QQ cuts []. Accordingly 
we conclude that we can identify and isolate the troublesome massless cut contribution to the 
two-point function. Indeed Eqs. ([[(]) and ([LI]) justify our previous Eq. (f|). This assertion 
might seem obvious but it is not : A QQ cut on the right-hand fermion loop in Fig. 0(a) does 
not provide, by itself, a pure real contribution. Only when both QQ cuts, on the left-hand and 
right-hand fermion loops of Fig. 0(a), are added we get ImllQQ. 

Let us go back then to Eq. (^). All the difficulty with the phenomenological application of 
the sum rules is now the fact that the contribution from the three-gluon cut is contained in 
both sides of the equality. Thus we propose an infrared safe definition of the moments by the 
trivial subtraction : 

^.^ crYi A ^ T-m Tl ( r*\ /v^ /"on A r^ Tm TT — lei 
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(13) 

Of course Eqs. flT2] ) and flT3| ) are meaningless unless we give a precise prescription about how to 
subtract the contribution of the massless cuts represented by ImII 3r Our previous discussion 
gives us the tool to proceed. Once the full 0(a^) IT (s) is calculated we can extract the imaginary 
part starting at s = (which should go with a 6(s) function) for any value of s. It is clear 
that the 6(s) and 8(s — AM 2 ) terms in the imaginary part of the vacuum polarization function 
correspond to three-gluon massless and to QQ cut graphs, respectively, and Im n 3g and Im IIqq 
are easy to distinguish, as Eq. (|TTD prevents the appearance of mixed 9{s) ■ 6{s — AM 2 ) terms. 
Therefore we identify Im Il3 5 and we now plug it in the dispersion integral of the right-hand 

4 The integration in Eq. ([ill) extends to the available three-gluon phase space. 
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Figure 2: Examples of perturbative non-heavy quark current correlators at 0(a 2 ) (a) and 0(a>l) 
(b) that contribute to the production of QQ states. 



side of Eq. (13) and perform such integration. Divergences contained in both this integral and 
M. n as q 2 — > will cancel with each other if the same infrared regularization is employed in the 
two quantities. The intuitive choice would be a low-energy cutoff s > 0, and Eq. (|13"D would 
be more precisely written as: 



M n = lim 

s ^0+ 



-7 (-4) n(g 2 



IV. 



dq 



-so 



1 r°° ds Imn 39 (s) 
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(14) 



where a vanishing therm in the So —>■ + limit has been omitted. 

The evaluation of the M. n moments at q 2 = < AM 2 made sense because, up to 0(a 2 ), 
this point is unphysical and the moments are well defined through an analytic continuation 
from the high-energy region. However note that the absorptive three-gluon contribution starts 
at q 2 = where perturbative QCD is unreliable. This introduces a further new difficulty in 
evaluating M. n moments at q 2 = 0, as we reach the physical non-perturbative region. Our 



definition of the moments, Ai n in Eq. (13), skips this problem by fully eliminating the massless 
terms and, therefore, the final heavy quark sum rule will only involve physics at q 2 > AM 2 . 

The method discussed in this Section could be extended to more general two-point correla- 
tors involving intermediate QQ states in their perturbative expansion. Indeed the correlator of 
two light quark vector currents has 0(a 2 ) contributions with an internal loop of heavy quarks 
(Fig. [^ (a)) and, similarly, the asymmetric correlator of a heavy and a light vector quark cur- 
rents is no longer vanishing at 0{a z s ) (Fig. ^ (b)). The absorptive part coming from the QQ 
cuts in the previous examples contribute to the phenomenological input o~(e + e~ — > QQ) in 
the usual sum rule analysis. Correspondingly they should be accounted for in the theoretical 
side. In short, the production of QQ states concerns not only the correlator of a couple of 
heavy quark currents and, for a more rigorous use of the sum rules method, this imbalance 
should be taken into account and properly fixed. A two-point function built from the sum of 
the electromagnetic currents associated to each quark flavour could be used in a generalized 



version of Eq. ([[]) : 

n£(p) = ^/rf 4 xe^^e (? e (? ,(0|T(g(x)7 M g(x))(7(0)7,g / (0))|0) , (15) 

where now q and q' stand for heavy or light quarks indistinctly, with electric charges e q and e q i, 
and q = q' is also allowed. As the different absorptive cuts contribute additively to Imll(p 2 ), 
the unwanted light quark and gluon qq, qqg, ggg, . . . cuts could be identified and, through the 
dispersive technique, subtracted from the full H(p 2 ) result. Consequently we are left with every 
possible QQ intermediate state arising from vector current production. Notwithstanding, the 
feasibility of this procedure from a technical point of view appears rather cumbersome and, at 
present, the experimental accuracy in the measurement of a(e + e~ — > QQ) cannot accommodate 
the corrections just discussed. 

4 Discussion on QCD sum rules applications 

The general rule given above is valid for all orders of perturbation theory, but it strongly 
relies in our ability to extract the massless absorptive part from the full result of IT(g 2 ) calculated 
at a definite order. Beyond 0(a 2 ) complete analytical results for the heavy quark correlator 
would be cumbersome and only numerical approaches may be at hand. In this sense, it would 
be convenient to have a method to calculate Im Hqq only based on Feynman graphs. We have 
already sketched such a method in the discussion following Eq. (|T0|) : we just need to sum up 
all the massless cut graphs to get Imn 3s , and then proceed with the dispersion integration that 
gives the associated dispersive part ||. For example, at 0(a%), the only massless absorptive 
part comes from the three-gluon cut in the diagram of Fig. |](a); let us call M.% g the amplitude 
producing three gluons from the heavy quark current at lowest order (i.e. through the quark 
triangle loop in Fig. [3]). The massless contribution to the absorptive part of the correlator is 
then: 

lmU 3g (s) = -j-JdR 3g M^ g -Mt gfl , (16) 

with the three-gluon phase space integral defined as 

1 1 IX 2 f s fS-ai 

dR 39 = ^tt^it: /„ ds ± /. ds ? ' ( 17 ) 
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in terms of the invariants s\ = (k\ + A^) 2 = {q — k^) 2 and S2 = (&2 + k^) 2 = (q — ki) 2 , and k{ 
being the momenta of the gluons. The real part would be obtained by integrating Eq. fllBT) : 



£o r- ds lmU 39 (s) = -so f- ds rs p: 
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which, in principle, could be performed also numerically. The nth-derivative of relation Qq ) 
respect to sq, in the limit So — > + , would give the infrared divergent contribution that should 
be subtracted from the full moments, as dictated by Eq. (|H|). 
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Figure 3: Feynman diagram for the production of three gluons at 0(a 



This evaluation solves the extraction procedure on the theoretical side. However the use 
of moments on heavy quark sum rules involves global quark-hadron duality which translates 
into the supposed equivalence between the theoretical evaluation of the moments and the phe- 
nomenological input. The implementation of the latter is not a trivial task. The total ex- 
perimental cross section a(e + e~ — ► hadrons) can be split into two disjoint quantities : the 
cross section for producing hadrons with Q-flavoured states, and the production of hadrons 
with no Q-flavoured components. If the experimental set up was accurate enough to classify 
events into one of these two clusters, the first class would be the required ingredient for the 
phenomenological part of the heavy quark sum rule. However this separation, implemented 
in the theoretical side within perturbative QCD, is rather involved. Up to 0(a 2 ) there has 
not been any doubt, in the literature, that contributions to this side arise wholly from QQ 
cuts in the heavy quark correlator Hqq- The physical picture behind this assertion relies in 
the assumption of factorization between hard and soft regions in the quark production process 
and subsequent hadronization. The hard region described with perturbative QCD entails the 
production of the pair of heavy quarks, and the soft part of the interaction is responsible for 
the observed final hadron content. Although possible, annihilation of the partonic state QQ 
due to the later interaction is very unlikely, as jets arising from the short distance interaction 
fly apart before long-distance effects become essential. Consequently, each jet hadronizes to 
a content of Q-flavoured states with unit probability. As local duality is implicitly invoked, 
this picture is assumed to hold at sufficient high energies; hence perturbative corrections to the 
hard part are successively included through the heavy quark currents correlator. 

The same physical picture applied to the three-gluon cut does not allow us to conclude 
whether this piece should be included or not in the theoretical side of the sum rules. Confine- 
ment tell us that this intermediate state hadronizes completely into hadrons with a content of 
light and/or heavy quarks indistinctly. It is conspicuous that if we could disentangle the heavy 
quark hadronization, 3g — > QQ, we should include only this piece into the sum rule. Then the 
singularity at q 2 = would disappear because heavy quarks are produced starting at q 2 = AM 2 . 
However there is no way to sort out light and heavy quark production off three gluons and, 



therefore, if we extract this contribution from the heavy quark sum rules we are introducing an 
incertitude in the procedure because we make sure that there is no light quark hadronization 
but we miss the heavy quark production. It is easy to see that the induced error is small, 
due to the fact that three gluons hadronize mostly to light hadrons. On one side, in the very 
high energy region and following perturbative QCD with Np = 4, we have only a 1/4 = 25 % 
probability of finding a specified pair of heavy quarks produced. And this is a generous upper 
limit because when we go down in energy, phase space restrictions severely reduce the counting 
of heavy quarks. Hence we estimate that excluding the three-gluon cut we introduce a tiny 
very few percent error in the sum rules procedure. 

5 Conclusions 

We have shown that rigorous and straightforward results of the general theory of singu- 
larities of perturbation theory amplitudes provide all-important tools to extract the unwanted 
0{a1) three-gluon massless cut pointed out by Groote and Pivovarov from the vector current 
correlator of heavy quarks. We conclude that the appropriate procedure to obtain information 
about the heavy quark parameters should make use of the infrared safe corrected moments, 
defined in Eq. (O), that now indeed satisfy the modified sum rule : 

M n = -Ids q? U , (19) 
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where the right-hand side can be extracted from the heavy quark production cross section 
a(e+e-^QQ). 

Finally we have pointed out that, starting already at 0{a 2 8 ), the use of sum rules associated 
to heavy vector current correlators shows and imbalance between the phenomenological input 
in the dispersion relation and the perturbative two-point function. This is due to the fact that 
the cross section of production of QQ heavy quarks is contained not only in the correlator of 
two heavy quark currents but in those involving light quarks too. We have indicated how to 
improve the application of sum rules by constructing a generalized correlator of vector currents 
in Eq. ( |i~5|) and, afterwards, extracting all the perturbative information not related with the 
production of heavy quarks, as we did in detail for the three-gluon cut. 
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